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Abstract 

The Total Ozone Mapping Spectrometer (TOMS) 
database is subjected to an analysis in terms of the 
Karhunen-Loeve (KL) empirical eigenfunctions. The 
concentration variance spectrum is transformed into a 
wavenumber spectrum, Edk). In terms of wavenumber 
Ec{k) is shown to be 0{k~^^^) in the inverse cascade 
regime, 0{k~^) in the enstrophy cascade regime with 
the spectral knee at the wavenumber of barotropic in- 
stability. The spectrum is related to known geophysical 
phenomena and shown to be consistent with physical di- 
mensional reasoning for the problem. The appropriate 
Reynolds number for the phenomena is Re « 10^". 

Atmospheric mixing is effected at horizontal scales 
which are large compared with the scale height 
(r^ 10km), which with inhibition of vertical motion by 
planetary rotation and stable stratification contribuljeg 
to the two dimensional picture of atmospheric activityE'EI 
In addition to its essential role in meteorology, current 
interest in mixing is enhanced by its role in regard to the 
behavior of the Antarctica ozone holm and to the lack 
of such an effect in the Arctic.Q 

Our investigation is based on satellite records of the 
earth's ozone fields. We have analyzed fifteen years of 
daily ozone fields of thCpTOMS (Total Ozone Mapping 
Spectrometer) database.El Due to technical and natural 
causes each daily record contains gaps in the form of 
missing pixels. A typical snapshot appears in Figure 
1. The dark regions represent areas of missing pixels 
caused partly by technical failure and in part due to 
polar night (measurements are based on reflected light). 
Each record is stitched together from sixteen separate 
records obtained from south-north synchronous orbits 
that are taken in a twenty four hour period from the 
satellite. Nimbus. 

The availability of such a large data set recommends 
statistical analyses, and we focus on spectral properties 
of the ozone field. Although ozone production (equa- 
torial regions) and depletioix. (polar regions) result from 
complex chemical reactionsQ these represent relatively 
weak sources and sinks and we follow common prac- 
tice and regard ozone as a passive scalar. The variance 
spectrum of a scalar contaminant in turbulent flows has 
been recently reviewed by Sreenivasan.Q Functional es- 



timates for the concentration variance spectrum (in ho- 
mogeneous isotropic turbulence) follow from dimensional 
arguments based on those first given by Kolmogorov 
leading to the famous Ek {k^ = Ke^^^k~^^^ energy spec- 
trum for the inertial rangew _ 

Both Obukhov and CorrsinEI show that an inertial 
range can exist, in particular the variance per wavenum- 
ber of concentration, c, denoted by Ec{k), has the form. 



E^ik) = Cxe-'^^k-^^\ 



(1) 



where C is a dimcnsionless constant, e is the usual tur- 
bulent energy transport rate, and x =< k(Vc)^ > is 
the appropriate dissipation rate. Thus the concentra- 
tion spectrum appears tied to the corresponding velocu 
ity spectrum. A review of experimental observationsQ 
shows departures from the universal form (Q), except 
perhaps at very high Reynolds numbec*. For rela- 
tively small diffusional effects Batchelonlj has shown 
that Ec — O {k~^) ■ This has been shown to hold un- 
der less restrictive hypothesis.0'0 Except for a recent 
simulationllj confirmation of this result has been elusive. 
Predictions of anomolous scalings have also appeared.t^ 

Arguments leading to the above spectra are unaltered 
when applied to two dimensional turbulence. However, 
the interpretation of the cascade of energy represented by 
the spectrum Ek {k) requires some additional remarks. 
Both energy E, and enstrophy = J (V A u) dx, are 
inviscid invariants in two dimensions. As a result of this 
Kraichnan and Batchelorllil have shown that the Kol- 
mogorov spectrum, Ek, represents an inverse cascade 
of energy from smaller to larger scales, and that there 
also exists a second cascade from small k to large k 
given by the enstrophy fl{k) — CoXo^^^ (with log cor- 
rection) and hence an energy spectrum E(k) cx k~^ , 
where C'p-is a dimcnsionless constant and xo = i'(V[j)^ 
(also seeEj). Suppcwt for E = 0{k^^) comes from many 
direct simulations .Il3 However, recent very large scale 
simulations show substantial divergence in the 0(fc-5/3) 
inverse cascade range.llZl Observational data from the at- 
mosphere is not definitive and although a power law en- 
ergy spectrum is indicated in the enstrophy f^nge the 
exponent appears to lie between ^2 and — 3.E3 In par- 
ticular, Schoeberl and BacmeisteiO suggest that the ex- 
ponent is —2 down to scales in the 10 kilometer range! 

A difficulty in interpreting these results for Ec already 
appears. Focusing on k large, it might be supposed in 
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analogy with three dimensions, that Ec (fc) = O (k^^) 
i.e., it should follow the energy spectrum. On the other 
hand, the vorticity (a scalar) formally satisfies the same 
convection equation as does a passive scalar, and from 
this one might suppose that E^, = O {k~^) , the Batchelor 
spectrum. As will be seen shortly, neither of these holds 
for Ec in the atmosphere. 

Other possible scalings for Ec have appeared in the lit- 
erature. For quasi two-dimensional turbulence Falkovich 
and Medijedevy find E = O (k''^^^) for laKe k. 
Saffman,E3 in considering the Burgers equation,Eil ob- 
served that its solutions are nearly piecewise discontin- 
uous which leads to Ec = O • Pierrehumbert us- 
ing concepts from chaotic mixing has obtained a variety 
of scaliwK for Ec from both mathematical and physical 
models.t3 

Satellite images (see Fig.l) are clearly inhomoge- 
neous and a transformation to wavenumber concepts 
is required- As will be seen the Karhunen-Loeve(KL) 
procedurecj is ideally suitedJor this purpose. In par- 
ticular, the snapshot methodcJ considerably reduces the 
needed computation effort. However, the presence p£ 
gappy data required modification of the methodology.ES 
This, as well as an eittensive analysis of the results, ap- 
pears in Manin et alEl and a mathematical treatment is 
also given elsewhere.tZl 

To connect the usual wavenumber spectrum with that 
obtained from the eaipirical eigenfunctions we recall 
an earlier discussionBa The concentration fluctuation of 
ozone is denoted by c(x, t) . For purposes of later di- 
mensional reasoning we write the dimensions of c as 
dim[c] = m/P, where m refers to molecules (of ozone) 
per area since the data are two dimensional. The mean 
variance in the homogeneous case is given by 



- 1 



c dx 



£c{-k)dk 



Ec{k)dk. 



(2) 



Thus dim[fc] = m'^l^^ and dim[i?c] = m^l^^ . To treat 
the inhomogeneous case corresponding to the data we 
consider the correlation 

i^c(x,y) -<c(x)c(y) > (3) 

which from KL can be written in spectral form, 

i^c = ^A„Vn(x)^„(y) (4) 

n 

where {"(An} are the eigenfunctions of the operator Kc- 
The total variance is given by 



< 



c^(x)dx >= TrK = 2^A„. 



(5) 



An eigenvalue A„ represents the average variance allo- 
cated to the projection of c onto '0„. The summation 
(||) is the natural generalization of (||) to the inhomoge- 
neous case. Each A„ represents the variance in a state, 
thus generalizing £'c(k) and has the same dimensions. 



dim [A] = dim[f ] = rn^l 



(6) 



In Figure 2 we display in doubly logarithmic form An 
versus index n. As is seen the variance spectrum falls, 
to good approximation, on two different power laws 



An 



n < 35, a, = .85 ± .035 
n > 50, Qfo = -1.56 ± .022 



(7) 



The error bounds appearing in (Q) and are based only on 
the least squares fit to the data, and not on the methods 
used in arriving at the spectrum which appears in Figure 
2. The region of the knee, 35 < n < 50, will be discussed 
below. 

In order to relate n to fc we observe that in the homo- 
geneous case modes carrying variances larger than those 
at k can be counted in number, N, as 



N oik^, whence k oc N^/'^. 



(8) 



Before continuing this reasoning it is instructive to 
verify the accuracy and content of these relations. For 
this purpose we employ the inverse relation between 
wavenumber and length scale. Thus (|^) implies that 
the length scale, L„, of the n**^ eigenfunction bears the 
following dependence on the index 



(9) 



An informal perusal of the eigenfunctions themselves 
supports this relationship between characteristic length 
and index. To quantify this we have-eomputed the cor- 
relation length of each eigenfunctionEa and the result is 
plotted in Figure 3. It is clear from this figure that (|) 
provides an excellent fit to the data in the two asymp- 
totic regimes. The region of the knee is the only anomaly 
and it appears as a plateau in the figure and corresponds 
to just one scale. 



27r/fc» 



« 4000 km. 



(10) 



It is generally stated in the geophysical literaturec^ 
that the baroclinic instability gives rise to a pattern, of 
wavenumber roughly seven. I.e., the unstable pattern is 
made up of approximately seven pairs of cyclonic/anti- 
cyclonic motions. With some indulgence on the part of 
the viewer, eigenfunction ■044 shown in Figure 4 appears 
to have this property. Roughly speaking, each of the 
eigenfunctions in the range 35 < n < 50 shows this spa- 
tial arrangement. To explain the plateau in Figure 3 we 
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suggest an analogy with the von Karman vortex trail. 
In that case a period seven disturbance rpauires four- 
teen independent modes for its description.Ej In view of 
the nature of the results this would appear to be a rea- 
sonable explanation. 

The energetics of the atmosphere has its origin in solar 
heating. However, dynamical activity introduces as 
the lengthscale at which mechanical energy is supplied, 
and is thus the significant lengthscale of the problem. 
If e is the energy transport rate which characterizes the 
inverse energy cascade, then 



Xo = k^e 



(U) 



characterizes the enstrophy cascade to higher wavcnum- 
bers. Using estimates for t he p hysical parameters of the 
atmosphere and L* from ( [lO|) we obtain the Reynolds 
number, R « 10^°. 

We now return to the implications of the power laws 
for A„, (0) to the wavenumber spectrum. In keeping 
with customary practice we consider the variance per 
wavenumber Ec{k) = k£c{k). It follows from (0) and (^) 
that 



Ec{k) cx 



/c-2/3, k<k^ 
k . k ^ k-^ 



(12) 



The more precise exponents are entered for suggestive 
reasons. The high wavenumber exponent lies slightly 
outside the errpij bound (However, it should be noted 
that KraichnarJlJ actually-estimated the enstrophy fall- 
off as l/kr'ln^/^k, see alsoLi^i With the exception of the 
Saffman-Burgers spectrum,E3 theoretical predictions lie 
outside the above range. 

In view of the relationship (|l2|), we now consider the 
consequences of elementary dimensional analysis. A 
straightforward argument yields. 



Equation (twe ) implies 

f(h/h \ ^ i '^iik/k*) ^^"^j fc/fc* -C 1 



(13) 



(14) 



where Ci and Co are dimensionless constants. In these 
terms 



Ecik) 



c^xe^'/^fc^'^"'/^ fc/A:, <1 ^^^^ 
CoXo ^^^xk*k~'^ fc/A;, > 1 



where the first form is appropriate for the inverse energy 
cascade and the second form is in a form appropriate for 
the enstrophy cascade. It should be noted that fc, (or 
i*) does not disappear under either asymptotic limit. 



In fact as the functional form ( p^ implies it would be 
impossible to eliminate this parameter entirely in both 
limits. 

The spectra obtained above covers a wide range. The 
range fc/fc* < 1 extends tOf-tlae 10,000 km wavelength 
limit imposed by dynamicsEil'EI and the fc/Zc* > 1 ex- 
tends down to wavelengths of the order of IQOkm, the 
resolution of the data. Smith and YakhotEa have re- 
cently suggested that an inverse cascade for E (fc) exists 
for a range of wavelengths-.greater than 10 km (but see 
Schoeberl and Bacmeistenla) . This is based on the asser- 
tion that cumulus cloud activity acts as an energy source. 
We see no evidence for this but do not regard this as a 
contradiction since the resolution of the satellite data is 
greater than 100 km. It is of course vexing-, that much 
the theory and numerical experiments (alsotj) discussed 
earlier do not agree with the observed satellite analy- 
sis which we present above. Only the Saffman-Burgers 
spectrum shows agreement. To test this further we have 
looked at the 'discontinuity' patterns of the data and 
find that Iv^l^ shows filamentous one dimensional pat- 
terns possible explanation for such strand-like pat- 
terns has been discussed recently. Both satellite obser- 
vations and computer simulations show the presence of 
tongues of stratospheric air extending from the tropics to 
mid-latitudes. These result from the breaking of Rossby 
waves at the edge of the polar vortices.a 
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